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Motivated by recent experimental work, we use first-principles density functional theory methods 
to conduct an extensive search for low enthalpy structures of CeCa under pressure. As well as 
a range of buckled structures, which are energetically competitive over an intermediate range of 
pressures, we show that the high pressure system (> 18 GPa) is unstable towards the formation of a 
novel class of layered structures, with the most stable compound involving carbon sheets containing 
five- and eight-membered rings. As well as discussing the energetics of the different classes of low 
enthalpy structures, we comment on the electronic structure of the high pressure compound and its 
implications for superconductivity. 
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Although the history of graphite intercalation com- 
pounds (GICs) dates back more than a century, they first 
became prominent in the early '60s when superconductiv- 
ity was discovered in some alkali metal GICs [l[ . Interest 
in these compounds has been reignited by the recent dis- 
covery that at least two compounds, CeCa and CeYb, 
have superconducting transition temperatures that, al- 
though modest by the standard of cuprate compounds, 
are an order of magnitude higher than those found previ- 
ously for GICs 0, 03] ■ Soon after this discovery, electronic 
structure calculations revealed that a three-dimensional 
band, originating from the intrinsic interlayer band of 
graphite and a metal ion band of the intercalant, plays a 
crucial role in facilitating superconductivity in GICs [J] 
(see also Ref. Subsequently, first-principles calcula- 
tions for CeCa and CeYb [a, 0] provided evidence that 
the out-of-plane phonons of the graphene sheets and the 
in-plane phonons of the metal atom couple to the elec- 
trons in this three-dimensional band and provide a basis 
for understanding the superconductivity within a con- 
ventional BCS framework. A brief review of the work to 
date is given in Ref. [§]. 

As well as superconductivity, evidence for a reversible 
pressure-driven phase transition has been reported in 
both C 6 Caand C 6 Yb 0,0. As the pressure is increased 
from ambient, the superconducting transition tempera- 
ture increases markedly, rising in CeCa from 11.5 K to 
15 K at 7 GPa above which it drops abruptly to 5 K. 
Similar behaviour is reported in CeYb with a transition 
at ca. 2.2 GPa. Since the transition is reversible, and cal- 
culations have shown softening of a phonon mode with 
increasing pressure, it was conjectured that a structural 
phase transition takes place ll|, |l2j . The reversible na- 
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ture of the observed transitions suggest that they involve 
small atomic displacements 22]. However, the various 
observations motivate an investigation of the wider land- 
scape of structures of GICs under pressure, including the 
possibility of large-scale rearrangements of the atoms. In 
this letter we report a study of possible structures of 
CeCa under pressure, finding that large-scale atomic re- 
arrangements are favoured at quite low pressures. We 
find that the behaviour of the CgYb shares a similar phe- 
nomenology. 

Before discussing the methodology, let us first sum- 
marise the main findings. The results below reveal that 
the low enthalpy structures of CgCa (and other GICs) 
can be broadly classified into different classes. While the 
planar structure of CgCa is stable at ambient pressure, 
structures involving a buckling of the graphene sheets 
becomes energetically competitive over a range of inter- 
mediate pressures. At a comparatively low pressure of ca. 
18 GPa, however, we find that the familiar honeycomb 
lattice structure of the graphene layers becomes unsta- 
ble towards a rearrangement involving large- and small- 
diameter carbon rings, with the former accommodating 
the metal ion intercalate. With such an arrangement, 
the volume may be reduced without greatly increasing 
the internal energy, which results in a low enthalpy. 

To arrive at these conclusions, we have have carried out 
an extensive search within the "space of possible struc- 
tures" of CeCa by relaxing a large number of random 
structures at a constant pressure of 15 GPa. (The gen- 
eral methodology parallels that described in Ref. [13j|.) 
This approach allows a search of the structure space with 
unbiased initial conditions. The calculations were per- 
formed within density functional theory jdl| using the 
castep package [3| , with a plane wave basis set and ul- 
trasoft pseudopotentials. To remain consistent with our 
previous work on the electronic structure of CeCa Q , we 
chose to use the local density approximation [161 ] for the 
exchange-correlation functional and used the pseudopo- 
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tentials that come with the CASTEP package, dividing the 
electrons into 2 core electrons and 4 valence electrons for 
C, 10 core and 10 valence for Ca, and 42 core and 28 va- 
lence for Yb. We tested the relative stability of the most 
stable structures with a generalised gradient approxima- 
tion 17J and found no significant changes. The initial 
conditions for the geometry optimisations consisted of 
a unit cell containing one formula unit (7 atoms) with 
random atomic positions and random lattice vectors, the 
latter bounded so that the unit cell volume was within 
a factor of two of the ambient value. For the structure- 
space search we used a medium quality of plane wave 
cut-off energy (270 eV) and Brillouin zone integration 
(0.07 A -1 ), so that the optimisations were fast, thus en- 
abling a larger number of runs to be carried out using the 
available computing resources (32 processor Xeon clus- 
ter). The search is believed to be nearly exhaustive over 
the range of 7-atom cells because the low enthalpy struc- 
tures came up repeatedly during the search. The Pmma 
structure with a 14-atom unit cell was built by hand us- 
ing the intuition gained from the results of the random 
search. 

As mentioned above, the search produced a number of 
novel structures which can, broadly, be grouped into two 
families. The first consists of the well-known rhombohe- 
dral or trigonal Rim structure which is the stable low- 
pressure phase of CgCa, and its lower symmetry variants 
in which the bonding topology of the graphene sheets 
remains intact, but the sheets are buckled (see Fig. [TJ. 
In these buckled configurations, the Ca ions occupy the 
troughs between the buckled sheets, which results in a 
small reduction in overall volume. The principal differ- 
ence between the Cm and the C2/m structures is that, 
in the former, the metal sheets are stacked in an aflj se ~ 
quence, just as in the Rim structure, whereas the latter 
has an aa stacking. Structures consisting of nano-porous 
carbon framework filled with Ca atoms were also found, 
but they were higher in enthalpy at these pressures. 

Despite the buckling of the sheets, the integrity of 
the hexagonal carbon ring structure in this family of 
compounds is maintained. In the second family of low 
enthalpy compounds, the carbon and Ca atoms main- 
tain a planar geometry, but the hexagonal ring struc- 
ture is replaced by a network of different sized rings. 
The Pmma structure has five-, six- and seven- membered 
rings, whereas the Cmmm structure has only five- and 
eight-membered rings. A useful way to think about these 
bonding topologies is that they result from one (Pmma) 
or two (Cmmm) Stone- Wales (SW) bond rotations [HI, 
starting from the regular hexagonal graphene network. 
The last two structures are reminiscent of those of a class 
of ternary boride compounds, whose prototype is YCrB4. 
To our knowledge, its structure, consisting of metal ion 
layers and boron sheets with five- and seven-membered, 
was first proposed by Kuz'ma [19j , Several other mem- 



Space group Lattice parameters 

(A, °) 

Cmmm a=9.07 6=3.66 c=3.54 
a=90 /3=90 7=90 



Atomic coordinates 

(fractional) 
CI 0.079 0.000 0.000 
C2 0.175 0.308 0.000 
Ca 0.500 0.000 0.500 



R3m a=4.24 6=4.24 c=12.28 C 0.000 0.333 0.167 
a=90 /3=90 7=120 Ca 0.000 0.000 0.000 



Pmma a=4.82 6=3.74 c=6.7 
Q =90 (3=90 7=90 



CI 0.250 0.500 0.117 
C2 0.514 0.500 0.210 
C3 0.750 0.500 0.091 
C4 0.596 0.500 0.417 
Ca 0.250 0.000 0.350 



C2/m a=7.30 6=4.23 c=4.15 CI 0.000 0.334 0.500 
a=90 /3=102 7=90 C2 0.171 0.168 0.532 
Ca 0.000 0.000 0.000 



Cm 



=7.32 
=90 



4.24 
=105 



=6.4 CI 0.246 0.166 0.513 
=90 C2 0.054 0.332 0.474 
C3 0.413 0.333 0.512 
Ca 0.531 0.500 1.000 



TABLE I: Details of CeCa structures at 15 GPa, calculated 
at the higher level of precision. 



bers of the class have since been identified 2fj, 21 1 . 

To further increase the accuracy of the procedure, we 
chose promising representative structures from the dif- 
ferent families and refined them using a higher plane 
wave cut-off energy (480 eV) and Brillouin Zone sam- 
pling (0.05 A -1 ). The converged lowest enthalpy struc- 
tures obtained are defined in Table [U Figure [2] shows the 
enthalpies of the new structures, referenced to the Rim 
structure, as a function of external pressure. The buckled 
Cm structure is slightly less stable than the Rim struc- 
ture at pressures up to about 18 GPa, but at higher pres- 
sures it becomes more favourable. In fact the Rim and 
buckled structures are always close in enthalpy, and their 
relative stability cannot be conclusively asserted because 
of inherent systematic errors in the local density approx- 
imation (commonly estimated as 0.05 eV/atom for total 
energy differences). The Pmma and Cmmm structures 
show a steep decline in relative enthalpy, indicating that 
they would be favoured at high pressure. Considering 
the geometry of these structures, it can be seen that the 
energy cost of the SW bond rotations is offset by a signif- 
icant reduction in volume (up to 20%), as the metal ions 
are accommodated within the larger rings. The enthalpy- 
pressure curve of Pmma has a lower slope than that of 
Cmmm because it involves one SW bond rotation rather 
than two, so that the rings in which the Ca atoms sit are 
smaller and the volume reduction is not as large. 

Figure[3]shows the band structure corresponding to the 
Cmmm structure. Comparing it to that of the "empty" 
carbon skeleton (after removal of the metal ions), it can 




FIG. 1: Low enthalpy structures of CeCa found during the 
structure space search. The top panel shows the familiar R3m 
structure, which is stable under ambient conditions. The Ca 
ions (large) are situated between the layers, occupying sites 
above the centres of the hexagonal rings in an eta stacking 
arrangement. The following two panels show lower symmetry 
variants, with C2/m and Cm structures, where in each case 
the graphene sheets are buckled (as depicted in the side views 
shown as insets) and the Ca ions are rearranged. While all 
of these buckled structures are energetically competitive, it 
is not possible to conclude from our calculations whether any 
become globally stable at intermediate pressures. The bottom 
panels show structures in which the hexagonal carbon rings of 
the graphene sheets are transformed, by a sequence of Stone- 
Wales bond rotations, into five-, six- and seven-membered 
(Pmma) or five- and eight-membered rings (Cmmm). At high 
enough pressure, the latter structure is expected to become 
globally stable. 
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FIG. 2: The enthalpy differences per CeCa formula unit of 
our new structures, referenced to the R3m structure. 



be seen that, upon intercalation, a new dispersive band 
becomes occupied in much the same way as in the parent 
R3m structure 0]. However, the density of states (DoS, 
shown in Figure[4]) at the Fermi level is significantly lower 
(on both the Ca and C atoms) in the Cmmm structure, 
suggesting that the conditions for superconductivity are 
likely to be less favourable. In addition, the general shape 
of the DoS looks much less like that of graphite as com- 
pared with the R3m structure. 

Taken together, these findings present a coherent pic- 
ture of the behaviour of C^Ca under pressure. At low 
pressures, the graphene sheets remain flat, leading to the 
familiar Rim structure. However, this parent structure 
is but one member of a larger family of compounds in 
which the integrity of the hexagonal rings is maintained, 
while the graphene sheets become buckled to accommo- 
date Ca atoms within the troughs. The troughs allow 
more space for the Ca atoms, so these structures are 
favoured over i?3m at higher pressures, although the gain 
is limited. At still higher pressures, another more drastic 
solution emerges in which the bonding topology of the 
carbon atoms is disrupted by SW bond rotations within 
the graphene layer. A single SW bond rotation results 
in five-, six- and seven-membered rings, as in the Pmma 
structure, and two SW bond rotations result in five- and 
eight-membered rings as in the Cmmm structure. The 
SW mechanism results in large-diameter carbon rings, 
within which the Ca atoms sit, giving a substantial vol- 
ume reduction, and these structures become stable at 
pressures above about 18 GPa. The accompanying phase 
transition is expected to be strongly first order, and the 
electronic structure of the large-diamctcr-ring structures 
is strongly modified from that of RZm. In summary, and 
recalling the experimental evidence for a phase transi- 
tion, we conclude that our simulations strongly predict 
a transition to the structure like Pmma or Cmmm at 
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FIG. 3: The band structure of the Cmmm structure (left), 
and the corresponding "empty" structure with the metal ions 
removed (right). The two panels have been aligned vertically 
using the bottom of the sigma bands (not shown), the respec- 
tive Fermi levels are indicated by dashed lines. Since the space 
group of this structure is almost hexagonal, we have used the 
corresponding notation for the points of high symmetry in 
the Brillouin Zone and indicated the small discrepancy by us- 
ing primed letters. Note the increase in the Fermi level upon 
intercalation, accompanied by the lowering and subsequent 
occupation of the dispersive interlayer band (marked by an 
arrow) . 



3 













Total 

-c 

Ca 




Cmmm \J ■ \ 




Total 

---C 
Ca 




R3m 


j. ,-- V f V\ 



-5 -4 -3 -2 -1 1 2 3 4 5 
eV 



FIG. 4: Density of states, resolved into contributions on the C 
and Ca atoms, corresponding to the Cmmm (top) and R3m 
(bottom) structures, calculated at 15 GPa. 

sufficiently high pressures, but due to the intrinsic er- 
rors associated with DFT, it is uncertain whether the 
observed transition is this one. It is possible that one 
of the buckled structures becomes the most stable one 
for an intermediate range of pressure, in which case we 
predict a second phase transition at higher pressure. 

Further calculations show that the bonding between 
the metal and carbon atoms is stronger than might be 
expected, and all the structures considered here are sta- 



ble against phase separation into diamond/graphite and 
pure fcc/bcc Ca. We have also studied the stability of 
above phases of CeYb under pressure, and again we find 
that structures with Yb atoms sitting within larger mem- 
bered rings are stable at pressures above about 18 GPa. 
We conjecture that the occurrence of such rings accom- 
modating the intercalate atoms might be a general fea- 
tures of graphite intercalation compounds under high 
pressures. 
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